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ABSTRACT 

Context. Recent progresses and constraints brought by helio and asteroseismology call for a better description of stellar interiors and 
an accurate description of rotation-driven mechanisms in stars. 

Aims. We present a detailed analysis of the main physical processes responsible for the transport of angular momentum and chemical 
species in the radiative regions of rotating stars. We focus on cases where meridional circulation and shear-induced turbulence only 
are included in the simulations (i.e., no internal gravity waves nor magnetic fields). We put special emphasis on analyzing the angular 
momentum transport loop and on identifying the contribution of each physical process involved. 

Methods. We develop a variety of diagnostic tools designed to help disentangle the role of the various transport mechanisms. Our 
analysis is based on a 2-D representation of the secular hydrodynamics, which is treated using expansions in spherical harmonics. 
By taking appropriate horizontal averages, the problem reduces to one dimension, making it implementable in a ID stellar evolution 
code while preserving the advective character of angular momentum transport. We present a full reconstruction of the meridional 
circulation and of the associated fluctuations of temperature and mean molecular weight along with diagnosis for the transport of 
angular momentum, heat and chemicals. In the present paper these tools are used to validate the analysis of two main sequence stellar 
models of 1.5 and 20 Mq for which the hydrodynamics has been previously extensively studied in the literature. 
Results. We obtain a clear visualization and a precise estimation of the different terms entering the angular momentum and heat 
transport equations in radiative zones of rotating stars. This enables us to corroborate the main results obtained over the past decade 
by Zahn, Maeder, and collaborators concerning the secular hydrodynamics of such objects. We focus on the meridional circulation 
driven by angular momentum losses and structural readjustements. We confirm quantitatively for the first time through detailed 
computations and separation of the various components that the advection of entropy by this circulation is very well balanced by the 
barotropic effects and the thermal relaxation during most of the main sequence evolution. This enables us to derive simplifications 
for the thermal relaxation on this phase. The meridional currents in turn advect heat and generate temperature fluctuations that induce 
diff'erential rotation through thermal wind thus closing the transport loop. We plan to make use of our refined diagnosis tools in 
forthcoming studies of secular (magneto- )hydrodynamics of stars at various evolutionary stages. 
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1. The impact of rotation on stellar evolution 

Rotation, and more precisely differential rotation, has a major 
impact on the internal dynamics of stars, in several ways. It in- 
duces large-scale circulations both in radiative and convective 
zones that advect simultaneously angular momentum, nuclides 
and magnetic fields (,Eddington||T925| [Wgt|[T925t [Sweet,, 1950( 
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When the star rotates differentially, various instabilities de- 
velop (secular and dynamical shear instabilities, baroclinic and 
multidiffusive instabilities) that generate hydrodynamical turbu- 
lence in radiative zones, in addition to these circulations. Just as 
in the terrestrial atmosphere and in laboratory experiments, this 
turbulence acts so as to reduce its cause, namely the gradients 
of angular velocity and of chemical composition. This explains 
why its effect may be described as a diffusion proces s ([Talon &| 
Zahn 1997, Maeder 2003] [Mathis et al.||2004t |Talon||20071 and 
references therein). 

Rotating stars have an equator that is cooler than the poles, 
which has an important effect on radiatively driven stellar winds 
and hence on the loss of mass and angular momentum ( [Maeder[ 



2 



Decressin et al.: Secular Hydrodynamical Processes in Rotating Main Sequence Stars 



|1999| l. If rotation is large enough, the stars can even reach the 
break-up limit, when centrifugal forces balances the gravity, and 
matter is ejected through an equatorial mechanical wind which 
create a circumstellar disk ( Maeder & Meynet|2000a[[Meynet et| 
|al.|2007t|^s55met al.,2008l [ 

Rotation also has a strong impact on the stellar magnetism. 
In radiative zones, rotation interacts with the magnetic field 
and is able to trigger magnetohydrodynamical instabilities, that 
could play a role in the transport of angular momentum and 
nucUdes ([Ch arbonnea u & MacGregOT] [1993] IGaraud"2002a 
Spruit|1999| 12002, Meno u et al. |2004 [ [Maeder & Meynet 2004{ 



Eggenberger et al. 2005 ; Braithwaite & Spruit 2005 ; Mathis & 



Zahn 2005, Braithwaite ,2006, Brun & Zahn 2006; Zahn et al. 
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Furthermore, internal gravity waves and gravito-inertial 
waves that are excited at the edge of the convection zones may 
also contribute to the transport of angular momentum. They 
propagate inside radiative zones and they extract or deposit an- 
gular momentum in the region where they are damped, thus 
modifying the angular velocity profile and the vertical distri- 
bution of chemicals (Press 1981; Schatzman||1993[ ITalon et al.l 
2002 , Talon & Charbonnel 2003, 2004 2005 2008 ; Charbonnel 
& Talon,2005,, Rogers & Glatzmaier,2005„Pantillon et al.,2007( 



Mathis et al.|2008r 

All these efi'ects significantly modify the evolution of rotat- 
ing stars. They affect their surface velocities and chemical abun- 
dances, and change their paths across the colour-magnitude dia- 
gram ( [Maeder & Meynet|200()bHMaeder|2"009l ). They also mod- 
ify the internal structure of stars in a way that we will soon be 
able to test thanks to asteroseismology over a broader range of 
evolutionary phases. 

In this paper we present a set of diagnosis tools adapted to the 
analysis of stellar evolution with rotation, and use them to com- 
pare the efficiency of different transport processes. The present 
study focuses on "type I rotational transport", where magnetic 
fields and waves are not accounted for, and where the angular 
momentum and the nuclides are both transported exclusively by 
meridional circulation and shear-induced turbulence. These di- 
agnosis and tools rely on a specific expansion of the equations 
for the transport of angular momentum, heat and chemicals that 
are briefly recalled in § 2. In § 3, we validate our method by 
applying our tools to two specific stellar models for which the 
hydrodynamics has akeady been extensively studied in the liter- 
ature. Conclusions and perspectives are presented in § 4. 



2. Modelling the secular processes 

2.1. Scale separation 

To simulate in full detail the dynamical processes in a star, would 
require to include length-scales and time-scales spanning many 
orders of magnitude. This is clearly not feasible, even with the 
most powerful computers. Either one chooses to describe what 
occurs on a dynamical time-scale, such as a convective turnover 
time, or one focuses on the long time evolution, as we shall here, 
where the typical time is either the Kelvin-Helmholtz time or 
that characterizing the dominant nuclear reactions. The same 
is true for the length-scales, at least in the vertical direction, 
where we shall take the resolution which represents adequately 
the steepest gradients that develop during the evolution. 

The situation is somewhat different in the horizontal direc- 
tion - i.e. in latitude, since we consider here only the axisym- 
metric case. Stellar radiative zones are stably stratified regions, 
and the buoyancy, which is the restoring force, acts to inhibit tur- 
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Fig. 1. Modelling strategy to study dynamical stellar evolu- 
tion. The diagram presents timescales of the typical physical 
processes as a function of the angular resolution necessary to 
properly describe these processes. The angular resolution is 
expressed in terms of the / index of the spherical harmonics 
Yi,m (0, (p). Inum - 600 indicates the maximum angular resolu- 
tion (in term of spherical harmonics nodes) presently achieved 
in global numerical simulations. 

bulent motions in the vertical direction. This leads to a strongly 
anisotropic turbulent transport that is more efficient in the hor- 
izontal direction (along isobars) than in the vertical one. As a 
result, the horizontal gradients of all scalar fields (temperature, 
angular velocity. . . ) are much smaller than their vertical gradi- 
ents, which allows their expansion in a few spherical harmonics. 
This scale separation, both in space and in time, is illustrated in 
Fig.g 

Considering the axisymmetric case, each scalar field is thus 
written as the sum of its horizontal average on an isobar, and its 
associated fluctuation, which is expanded the Legendre polyno- 
mials basis Pf(cos Q), up to some maximum order /",/,. For sim- 
plicity's sake, we shall present here the results keeping a single 
scale 1-2, that is dominant. Let us note that there are circum- 
stances where we have to include many more components, for 
instance when dealing with magnetic fields. 

The unresolved scales intervene in the turbulent transport, 
for which a prescription is applied that is derived, whenever pos- 
sible, from laboratory experiments or numerical simulations, and 
if not through phenomenological considerations. 

2.2. Linearization and expansion in spherical harmonics 

Let us briefly recall how these expansions are performed in prac- 
tice; for more details we refer the reader to Mathis & Zahn[ 
( 2004| l and [Maeder (2009 1. We begin with the macroscopic ve- 
locity field, which is split in 3 components: 



V ^ rsin6n(r,6)e^+ re, +t4M(r,0) 



(1) 



Term 1 represents the azimuthal velocity field associated with 
the differential rotation. Term 2 corresponds to the radial 
Lagrangian velocity due to the structural readjustments of the 
star during its evolution. Term 3 is the meridional circulation 
velocity field, "e^ where k - {r,0,(p} are the unit vectors in the 
radial, latitudinal and azimuthal directions. 
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In radiative zones, because of the strong stratification, the 
turbulent transport is highly anisotropic, which leads to a "shel- 
lular" rotation profile (Zahn 1992 ). The angular velocity thus 
varies little on an isobar, allowing to expand it as 



2.3.1. Transport of angular momentum and thermal-wind 

We first consider the azimuthal projection of the momentum 
equation. By averaging it over an isobar, we obtain the following 
advection-diffusion equation for the mean angular momentum : 



Q(r, 0) = Q(r) + Q2(r) [i^^Ccos 0) - 1/5] , (2) - 1 / ^ 1 /-^4^^ ^ ^ 1 (py,,r%m . (11) 

df ^ ' Sr-^ ^ / \ I 



where 



Q(r) 



j^^sin^ede 



(3) 



As usual, the meridional velocity field is projected on spher- 
ical functions, of which we show here only that of lowest order: 

d/*"* (cos B) 

'Um = Vi (r) Pi (cos 0)% + V2 (r) -\ ' eg , (4) 

w 

where f 2 is the second order Legendre polynomial. The anelas- 
tic approximation is adopted, filtering out the acoustic waves, 
which is amply justified for this slow meridional circulation. 
Therefore, the continuity equation reduces to V ■ {pHu) - 0, 
which leads to the following relation between the latitudinal (V2) 
and the radial (Ui) amplitude functions: 



(5) 



Moreover, since JiUm is divergenceless, it may be expressed us- 
ing a stream function (^m('", ^)) as in Spiegel & Zahn (1992), 
where 

P'^M^-TTTTT^V^x'e^ (6) 



r sin0 



that leads to 



P'Z/m;i- = — 7 ^. rP e^yv and p'Z/M;e = — r— t^^^m (7) 
r'- sin Q 



r sinfl 



with 



^jyi - -pr^U2 (cos^ 9 - cos 6^ . 



(8) 



P'Um, and thus Um, are then tangent to the iso-^ lines since 
1/MV^ = 0(c./Eq.|6). 

Finally, all other variables are expanded as shown here for 
the temperature and the mean molecular weight: 

T (r, e)^f (r)+6T (r, 6) with 6T (r, 0) = ["¥2 (r) f] P2 (cos 9) 
and 

H (r, 0)^J1 {r)+5ii (r, 0) with dfi (r, 0) = [A2 (r) /7] P2 (cos 0) ; 

(10) 

T and p are the horizontal averages, ST and dfx being their fluc- 
tuations and *P2 and A2 their relative fluctuations on the isobar. 



2.3. The transport equations 

Using these notations, the transport equations for angular mo- 
mentum and chemicals may thus be expanded in the following 
way. 



Term 1 represents the transport of angular momentum by the 
meridional circulation; note that the advective character of that 
transport is preserved. The diffusion term 2 is associated with 
the action of the shear-induced turbulence, where v,, is the tur- 
bulent viscosity in the vertical direction (see below Eq.[28]l. The 
Lagrangian time derivative d/dt is defined as d/dt - dt + rd,-. 

This equation shows the relation between the meridional cir- 
culation and the transport of angular momentum. In the asymp- 
totic regime, the left-hand side term is zero and the transport 
of angular momentum by meridional circulation is exactly bal- 
anced by that through shear turbulence. In the limit of vanishing 
turbulent viscosity, the rotation profile would adjust so that no 
meridional currents appear (Busse 1982V In reality because of 
the loss of angular momentum by the wind, and/or its redistribu- 
tion by structural changes, the left-hand side term is non zero. 

By integrating Eq.([TT]i over a surface of radius r, it can be 
recast in an equation for the fluxes: 



where 



r(m) = -Ftot = -Fmc (r) - Fs (r) , 



^^MC (r) = -ipr^nj/z 



is the flux transported by the meridional circulation, and 



(12) 



(13) 



(14) 



is that carried by the vertical shear induced turbulence. The term 
T{m) represents the loss or gain of angular momentum inside 
the isobar enclosing the mass OT(r); its derivation is given in 
Appendix [a] 



r(m) 



1 d 

47T dt 



Jo ' 



r'^Qdm' 



(15) 



From Eq. ^V2\ it is also possible to extract the radial component 
of the meridional circulation velocity: 



t/2 = f/r + Uv 



pr'^D. 



\r(m)- 



(16) 



Finally, by taking the curl of the momentum equation and again 
filtering out the fast timescales, we obtain the so-called thermal 
wind equation: 
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■ (^A2 - 5^2) 



(17) 



where 'g is the horizontal average of the effective gravit}{^ The 
coefficients 5 - - (dlnp/dlnT)p^ and ^ = {dlnp/dlnfijpj are 
introduced by the generalised equation of state (Kippenhahn & 
|Weigert| 19901 see |Maeder & Zahn|1998| or |Mathis & Zahn,2004 
more details). The right-hand side term is the barocUnic torque 



' In a rotating star, the centrifugal force diminishes the local gravity, 
and the effective gravity g is then given by horizontally averaging the 
hydrostatic equilibrium : g = VO -1- ^QrV (r sin d)-, where O is the 
gravitational potential. 
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2.3.2. Thermal relaxation 

Likewise, by expanding the energy conservation equation over 
spherical functions, we estabhsh the following equation for the 
temperature perturbation represented here by (Eq. 101 in 
IMathis & Zahn|2004| see |Zahn|1992l i: 



pTC, 



■ df 



- L 



_L _L 

- p—'2,T'n + p — '2,n-G 



M 



M 



pTO,S 



pTu,a,s 



pe 



(18) 

where Nt, the buoyancy frequency linked with the entropy strat- 
ification is given by: Nj = (gd/Hp) (Vad - V) with the usual no- 
tations for the temperature gradients: V = dlnT /dlnP. This is 
again an advection/diffusion equation, where the advective term 
plays the role of a heat source or sink, since the meridional cir- 
culation is caused by the transport of angular momentum, with 
the other terms describing essentially the thermal relaxation due 
to radiative damping. V ■ Fh is the entropy flux carried by the 
horizontal turbulence. 

For diagnostic purpose (and also for historical reasons; see 
the discussion in Zahn 1992 after Eq. (3.32)) we split the 



difiiision term in two pieces, which expressions are given in 
Appendix [B| 

- 7"2,s corresponds to the divergence of the radiative flux as- 
sociated to the deformation of the isobar induced by the per- 
turbing force, which is here the centrifugal acceleration. This 
term subsists in the case of solid body rotation; for this rea- 
son we called it the barotropic term following Zahn ( 1992| l. 

- Ta.Th contains the highest derivatives of the difiiision opera- 
tor and according to the thermal wind equation ( [T7] l it van- 
ishes in an homogeneous star when the rotation is uniform. It 
is referred to as the thermic or baroclinic term ( |Zahn|1992] l. 

The last term 7~2.n-g is associated with the nuclear reactions and 
radial structural readjustments during its evolution and usually 
plays a negligible role. 

Eq. fW\ can finally be written in the following form 



where 



S - TAdv 


+ T'B + Tjh + Tn-g, 


(19) 




- d^Pz 
-pTCp — , 


(20) 


'7~Adv - 


- Nl 
--pTC„f^U2, 


(21) 


t'B 




(22) 






(23) 




= P^nN-G. 


(24) 



This will be used in §3.4. in our diagnosis tools. 



2.3.3. Transport of nuclides 

The expansion of the transport equation for the nuclides on an 
isobar leads to an equation for the evolution of the mass fraction 
of each considered chemical (see also |Meynet £ Maeder]200()) : 



[ df ; 



d 



{Anr^pf{D, + D,s) 



dXi 



(25) 



where dMr = Anpr^dr and Z),, is the vertical component of the 
turbulent diffusivity (see Eq. [28] below). The strong horizontal 
turbulence leads to the erosion of the advective transport, which 
can then be described as a diffusive process (Chab oyer & Zahti] 
|1992| ) with the following effective diffusion coefficient 



eff 



30Dh 



(26) 



where D/, the horizontal component of the turbulent diffusivity 
(see Eqs. 29 and 30 1. 



The second term on the right-hand side of Eq. ( 25 1 coiTe- 
sponds to the temporal mass fraction evolution of the ("'"^nuclide 
due to nuclear burning. 



Equation ( 25 1 is complemented by an equation for the time 
evolution of the relative fiuctuation of the mean molecular 
weight, expressed here in terms of A2 : 



dA2 dlnu A^;^ 6 
at at gip r'- 



(27) 



where A^, the chemical part of the Brunt-Vaisala frequency, is 



given by A^ = {giplHp)V^ with = d\nJild]nP. 



2.3.4. Turbulence modelling 

The details of turbulence modelling have already been exten- 
sively discussed in previous papers ( [Talon et al.|1997[ |Palacios| 
et al.|200^ , so that we will just recall here the expressions we 
have chosen. 

For the vertical turbulent diffusion coefficients (D,, ^ Vy) we 
use the expression by Talon & Zahn ( 1997[ ): 



Nll{KT+Dh)+NllDh 



(28) 



where TJi'c = 1/6 is the adopted value for the critical Richardson 
number and Kj the thermal diffusivity. 

For the horizontal turbulent viscosity we have used two dif- 
ferent prescriptions: 



Vh 



■ ^[Crn\2V2 - at/21] where a = ^ 51n(rQ) 



d\nr 



after |Mathis et al.| ( [2004l ), and 

r 



Vh = —\2V2-aU2\ 



(29) 



(30) 



nucl 



after Zahn ( 1992 1. C - 1.6 x 10"*, Cr is a parameter of order 
unity, and V2 is given by Eq. (j5]). As for the vertical turbulent 
diffusion coefficient, we assume that Dh - v/,. 

We stress that the advective character of the transport of an- 
gular momentum by the meridional circulation makes the inter- 
pretation of the whole hydrodynamics more complex than when 
the diffusive approximation is used (e.g., End al & Sofia|[l978t 
Pinsonneault et al.|,1989 Heger et al.||2000) . This is one of the 
reasons that incited us to develop the tools that we describe in 
the following section. 



3. Numerical simulation of secular transport 

In this section, we apply the formalism described in § |2] to the 
evolution of rotating stars with initial masses of 20 and 1.5 Mq 
and metallicity close to solar value (Z = 0.017 and 0.02 for the 
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30000 20000 7000 6800 6600 6400 

Teff (K) Teff (K) 

Fig. 2. Evolutionary path in the HR diagram of a 20 Mq (left panel) and 1.5 Mg (right panel) star with solar composition. The solid 
line represents the main sequence. The dots on each plot indicate specific locations at which a detailed analysis of the structure is 
presented. 




0.5 1 1.5 

r (Re) 

Fig. 3. Angular velocity profile inside the 20 Mq {left panel) and 1 .5 Mq (right panel) model as a function of the (expanding) radius, 
at three different evolutionary stages on the main sequence, as indicated by the central hydrogen content. The short vertical lines 
indicate the location of the convection zones. In the 1.5 Mq model, at the end of the main sequence (Xc = 0), the convective core 
has disappeared. 



20 and the 1.5 M© respectively). In the Hertzsprung-Russell di- 
agram shown in Fig. [2] the black dots indicate the evolutionary 
points where a snapshot on the internal structure will be taken 
and analysed in details. These points are associated with central 
hydrogen mass fractions - 0.675, 0.32, 0.00 as indicated on 
the tracks. 

A number of results concerning the hydrodynamics of such 
rotating models have already been publish ed (Talon 1997 Talon 



& Charbonnel 1998, Meynet & Maeder||2000 
2003 ,Mathis et al.,,2004, ,Palacios et al.||2006[ 



Palacios et al. 
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2007 I, and we will rely on this experience to prove the rele- 
vance of the diagnosis tools that we have developed. The models 
presented here were computed with the stellar evolution code 
STAREVOL V2.90, and the reade r is referred to |Siess et al.| 
( |2000| i, [Palacios et al] ([2006} and [Siessl (|2006} for a detailed 



description of the input physics. We shall just recall the main 
characteristics and parameters used for the modeling. 



The reference solar composition we use follows Grevesse et 
[aLl ( |1996) . The standard mixing length theory is used to model 
the temperature gradient in the convection zones and the param- 
eter ffMLT - 1-75. The atmosphere is treated in the gray approx- 
imation, and integrated up to an optical depth t ^ 5 x 10"^. 

In both simulations mass loss is included from the Zero Age 
Main Sequence on. We use the |Reimers| ([T975 ) prescription for 
the low-mass star, with a parameter t/r = 0.5 (see also Palacios et 
al.|2006[ ) and the prescriptions by| Vink et al.| ( |2000) l for the mas- 
sive star Rotation effects on mass loss are accounted for follow- 
ing Alaeder_&_Aleynet( 2001). Angular momentum losses asso- 
ciated with mass loss are also accounted for, but their anisotropy 
( |Maeder|2002| l is not included in the computations. 
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Concerning the transport of angular momentum and nu- 
clides, we use Eq. ( [29] l as the prescription for D/, for the 20 Mq 
model. For the 1 .5 Mq model, we use Eq. ( 30 1 as in Palacios et al. 
( |2003| l who showed that it yields to a very good agreement with 
the surface velocity and the light nuclides (Li, Be) surface abun- 
dances that are observed in open clusters such as the Hyades. 
The expression for the vertical turbulent diffusion coefficient is 
given by Eq. (28 i in both models. For the transport of nuclides, 
we do not account for atomic diffusion. 

In our framework, the evolution of the angular velocity pro- 
file and of the meridional circulation in the radiative zones is 
governed by a system of five first order differential equations. 
These equations are obtained by splitting Eq. (Ill into first order 
equations that are complemented by Eq. ( 17 and Eq. (p7|. We 
use a Newton-Raphson relaxation method ( Heyney et al.| 1964) 1 
to solve these equations. In convection zones we assume solid- 
body rotation (0(r) = est) and do not solve the aforementioned 
system. The boundary conditions at the radiative zone limits are 
the same as those described in |Palacios et al.| ( [2003| l. 

We choose initial velocities on the ZAMS that reflect an av- 
erage velocity for main sequence B and F stars (Fukada 1982, 
Abt et al.|2002H Gaige 1993). The adopted initial velocity on the 
ZAMS are uzams = 300 km s"' and i/zams = 100 km s"' for the 
20 Mo and the 1 .5 Mq respectively. 

The evolution of the surface velocity can be affected by the 
magnetic torques exerted at the surface of the star during its early 
evolution ( Schatzman] 19621. This magnetic braking is gener- 
ally associated with the presence of a convective envelope ; it 
will thus only be applied to the low-mass star model, follow- 
ing | Kawaler's ( 1988 1 formalism as described in Palacios et al. 
( |2003) . In the case of massive stars, it is only very recently that 
magnetic fields have been detected in normal O and B stars, and 
they appear to be very weak or nonexistent ( Bouret et al.|2008 
[Schnerr et al.||20()8| ). Moreover, O and B stars with measured 
V sin ; appear to be fast rotators. For these reasons, we have de- 
cided not to account for magnetic braking in the 20 Mq model. 

In the following sections, we present new diagnosis tools to 
better understand and quantify the processes responsible for the 
building up of differential rotation in stellar interiors. 

3.1. Rotational profiles and differential rotation 
3.1.1. 20 Mo model 

The 20 Mq model does not undergo magnetic braking and thus 
remains a fast rotator during its main sequence (MS) evolution, 
with i;surf - 250 km s~'. Similarly to what was obtained by 
Meynet & Maeder ( |2000| l, differential rotation rapidly develops 
in the interior, as firstly meridional flow transports angular mo- 
mentum inwards. After a transient phase that lasts ~ 1 Myr fol- 
lowing the arrival on the ZAMS, the morphology of the rotation 
profile is not significantly affected until the star reaches the turn- 
off, as shown in Fig. [3] (left panel). 

Above the radial coordinate 5 Rq (resp. 6 Rq and 10 Rq) for 
Xc = 0.675 (resp. = 0.32 and X^ - 0.0) the angular velocity 
gradient becomes positive due to a very efficient transport of an- 
gular momentum by the meridional circulation in a region where 
the density is low (see § 3.2 1. This inversion of the angular veloc- 
ity gradient should be counteracted by the shear turbulence that 
tends to smooth out this profile. However, the shear can only 
operate if the flow is turbulent, i.e. if the Reynolds condition 
vv > vRcc is satisfied where Rcc - 1 is the critical Reynolds 
number and vy and v the shear and microscopic viscosities, re- 
spectively. Below the convective envelope (where - 0), 



© 




Fig. 4. Meridional circulation currents in the 20 Mq model when 
= 0.32. Solid thin (red) lines indicate counterclockwise cir- 
culation {U2 < 0), which carries angular momentum outwards. 
Hatched regions correspond to convection zones. 




Fig. 5. Decomposition of the meridional circulation based on the 
equation of angular momentum evolution in the 20 Mq star, 
when = 0.32. Hatched regions correspond to convection 
zones. Ur (short-dashed red line, see Eq. 16 1, U2 (meridional 
circulation, in long-dashed black line), and Uv (turbulent trans- 
port, in dotted blue line) are presented on a logarithmic scale. 
The thin part of the curves means that the plotted quantity is 
positive and the bold parts that it is negative. 



this condition is not fulfilled, and the shear is unable to connect 
the region with positive angular velocity gradient with the one 
where S^Q - (see the dip in the profile of Dy in Fig. 



15 1. This 



explains the persistence of this positive gradient. Computations 
done without taking into account the Reynolds criterion show 
that this positive gradient is reduced (by a factor 3 to 5) but is 
still present, although at a lower level. It should be noted that 
this inversion has already been obtained by Talon et al. ( I997| l 
in the asymptotic regime for a 9 Mq with an equatorial angular 
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velocity of 425 km s ' . This feature can be linked to the Gratton- 
Opik term in the meridional circulation (see § 3.4.1). 

The decrease in the angular velocity during the MS evolution 
is due to two factors: (1) mass and angular momentum losses by 
the stellar winds and (2) structural changes that lead to an overall 
expansion of the star It should be noted however that up to the 
point where - 0.1, the contrast in Q. between the center and 
the surface is small, not exceeding 20% . 

3.1.2. 1.5 Mo model 

Contrary to the 20 M© star, the 1 .5 Mq star undergoes magnetic 
braking at its arrival on the ZAMS, resulting in a strong extrac- 
tion of angular momentum, as shown in Talon & Charbonnel, 
( |T998| l and |Palacios e"taL] ( p003]|2006| . During the first 800 Myr 
the surface equatorial velocity drops from 100 km s ' to 78 km 
s in agreement with observations in open clusters like the 
Hyades (i.e., Gaige 1993). Beyond this point, the braking be- 
comes less efficient and the surface rotational velocity contin- 
ues to decrease slowly until the star reaches the turn-off with 
I'surf - 53 km s~' at age 2.35 Gyrs. 

The right panel of Fig. [3]presents the evolution of the angu- 
lar velocity profile on the MS. When = 0.675, i.e. when the 
star is ~ 200 Myr old, the core already rotates 1.5 times faster 
than the surface. While the angular velocity in the convective 
core slowly increases, the extraction of angular momentum at the 
surface slows down the surface layers and differential rotation 
increases during the MS. When the convective core disappears 
at central H-exhaustion, the contraction of the inner shells pro- 
duces a much steeper O-gradient (dashed profile in right panel, 

Fig.|3]l with Qcenter = 14 X Qjurf. 



3.2. Meridional circulation 

In this section, we isolate the processes which are at the origin 
of meridional circulation. 



3.2.1. 20 Mq model 

Figure |4] shows the meridional flows in the radiative interior 
of the 20 Mq when the central hydrogen mass fraction reaches 
Xc = 0.32. They are represented in terms of the stream function 
(r, 0) (see Eq.|8]l. The circulation consists of a single counter- 
clockwise loop, by which matter is transported inward along the 
rotational axis, and is conveyed outward in the equatorial plane. 
These flows thus extract angular momentum from the deep inte- 
rior to the surface. Matter is conserved in the meridional circula- 
tion currents and p H'J/mII = constant along the flux lines (shown 
in Fig. [4]). The abrupt jump in density occurring near the surface 
produces an increase in the meridional circulation velocity as 
illustrated by Fig. [5] This in turn force s the inversion of the an- 
gular velocity gradient mentioned in § 3.1 via the Gratton-Opik 

in Appendix Bl). The velocities 



—2 



B.l 



term (oc -Q /2nGp, see Eq. 
are maximum in the low density regions (i.e. below the surface) 
where they reach W'tluif, 0)\\ ~ 1.6 m s"'. 

After a short adjustment phase, which depends on the initial 
conditions that have been imposed on the ZAMS (uniform rota- 
tion), the meridional circulation settles into a counterclockwise 
regime, as was noticed by Mey net & Maeder] ( |2000 ). This flow 
thus carries angular momentum outwards, in order to compen- 
sate the structural changes, i.e. the inflation of the envelope and 
the contraction of the core. 



Figure |5] shows the a posteriori reconstruction of the merid- 
ional circulation based on the integration of the equation for the 
evolution of angular momentum (Eq. 16 1, at - 0.32. Ur rep- 
resents the angular momentum extraction and Uy that due to 
the shear. In the inner part of the radiative zone, below 7 Rq, 
the meridional circulation is powered by the extraction of angu- 
lar momentum, as \Ur\ » \Uy\. Near the surface, Uy becomes 
negative and of the same order as Ur- This inversion of Uy is 
correlated with the sign change of the angular velocity gradient 
since Uy = -5v,,5,. InQ. Ur may be considered as a measure 
of the departure from a stationary regime, where angular mo- 
mentum transport by shear (Uy) and meridional circulation (U2) 
compensate each other. This balance is barely achieved in the 
uppermost part of the star (above 8 Rq), but the bulk of the ra- 
diative interior is not in a stationary state, as already pointed out 
in [Talon et aL](|1997) and Meynet & Maeder (2000) on the basis 



of an asymptotic analysis. Here, we confirm this result but also 
track down the importance of the shear versus angular momen- 
tum transport in shaping the meridional circulation. 



3.2.2. 1.5 Mo model 

The wind-driven meridional flows are depicted in Fig. |6] for the 
low-mass star, at the three evolutionary stages indicated in Fig.|2 
Compared to the 20 Mq model, the topology of the meridiona 
circulation is more complex. In panels 1 and 3, the circulation 
presents two loops: a small clockwise loop in the central regions 
and a large counterclockwise loop connecting the interior to the 
surface. As in the massive star model, the counterclockwise loop 
carries angular momentum outward and is driven by the strong 
extraction of angular momentum (here largely due to the mag- 
netic braking). The clockwise loop corresponds to matter flow- 
ing from the equator to the pole, which results in deposition of 
angular momentum in this region. Such loops appear only in the 
center and do not persist over the whole evolution on the main 
sequence. 

At the beginning of the MS evolution, the meridional cir- 
culation velocity is maximal in a thin layer below the surface 
as a result of the applied magnetic torque. The amplitude of the 
meridional circulation velocity \\U2\\ does not exceed 0.5 m s"' at 
its maximum and is of the order of 10 in the bulk of the radia- 
tive zone. As the star spins down, the braking becomes less effi- 
cient and the meridional circulation weakens by 3 orders of mag- 
nitude. The large difference in meridional velocity between the 
low-mass and the massive star is due to the L/M factor (Eq. [T8]l, 
which increases with stellar mass, and to the higher initial angu- 
lar velocity of the 20 Mq star 

Figure [T] is similar to Fig. |5]but for the 3 models selected 
during the MS evolution of the 1.5 Mq star As for the 20 Mq 
model, the meridional circulation is driven by the local loss of 
angular momentum (Ur term which matches the Uo profile), 
the contribution of the vertical shear (Uy) being always much 
smaller, except in the central regions where it can be of the same 
order as Ur. This figure confirms that the inward transport of 
angular momentum (clockwise meridional circulation loops) is 
associated with a positive value of Ur, corresponding to a local 
gain of angular momentum (see Eq.[T6|. In most of the radiative 
zone, Uy <sc U2, and contrarily to what was found in the massive 
star model, the circulation in the upper part of the radiative zone 
never approaches the stationary regime. As shown by |Talonetal.| 
( 1997 ), the extraction of angular momentum by magnetic brak- 
ing together with the very small meridional velocities prevent the 
inversion of the angular velocity gradient in the regions of low 
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Fig. 6. Meridional circulation flow in the 1.5 Mq model when = 0.675 (left panel), 0.32 (middle panel) and 0.0 (right panel). 
Solid thin (red) lines indicate counterclockwise circulation (U2 < 0) and solid bold (blue) lines clockwise circulation (U2 > 0). In 
this model, the outer cell is turning counterclockwise allowing the extraction of angular momentum by the wind. Hatched regions 
indicate convection zones. 
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Fig. 7. Decomposition of the vertical component of the meridional circulation U2 according to Eq. ( 16 1 for the 1.5 Mq model at the 
three different epochs indicated on top of each panel. The term due to angular momentum transport. Up, is shown in dashed red line, 
the one due to vertical shear, Uy, is in dotted blue line. The value of U2 obtained in the simulation is given in long-dashed black 
line. Thin and thick lines represent respectively positive and negative values. 



density as was the case in the massive star model. Consequently, 
Uy oc -d,Q. remains always positive. Besides, during the main 
sequence the differential rotation remains too small to ensure an 
efficient transport of angular momentum by the shear turbulence. 

3.3. Angular momentum transport 



We already compared in § |3.2| the respective contributions of 
meridional circulation and turbulence in the transport of angu- 
lar momentum, in agreement with results from the literature. 
We shall now examine the problem from a different perspective, 
based on the fluxes carried by the two processes (Eq. 12 1. 



3.3.1. 20 Mq model 

Figure |8] shows the two components of the angular momen- 
tum flux associated with the meridional circulation and shear 
In this model the meridional circulation presents only one coun- 
terclockwise loop so the transport occurs in the same direction 



in the whole radiative zone, and carries angular momentum from 
the core to the surface (Fmc > 0). Below r - 6Rq the turbulent 
shear works in the same direction (F^ > 0), but its contribution is 
negligible compared to that of the meridional circulation. Above 
6 Rq, the angular velocity gradient changes sign and shear takes 
over the advection of angular momentum. 

and Fmc are of opposite sign below the convective en- 
velope during the MS evolution, and grow in amplitude due to 
the radial expansion of the outer low-density layers (see Eqs.p3] 



and 14 1. This configuration tends to impose a stationary circula- 



tion in this region. On the other hand, these fluxes have the same 
sign deeper in the interior, and act jointly to compensate for the 
structural changes of the star by redistributing the total angular 
momentum. 



3.3.2. 1.5 Mq model 

In the low-mass model, the meridional circulation dominates at 
all times the angular momentum transport and Fmc - ^Xot on 



Decressin et al.: Secular Hydrodynamical Processes in Rotating Main Sequence Stars 



X 

o 



= 0.675 



0.32 



0.0 




0.5 1 0.5 1 1.5 0.5 1 1.5 2 

r (Ro) r (R^) r (R^) 

Fig. 9. Same as Fig.jijfor the 1.5 Mq model at the three selected evolutionary stages indicated by X^. In this graph, the logarithm of 
the absolute value of the flux is plotted and rescaled by R^. Thin and bold lines represent positive and negative fluxes respectively. 
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Fig. 8. The total flux of angular momentum (solid black line) de- 
composed into its meridional circulation F^c ('") = -\pr''^^U2, 

(dashed red line) and turbulence Fs (r) = -pVyr^drO. compo- 
nents (dotted blue line) for the 20 Mq model when - 0.32. 
The fluxes are rescaled by R'^. 



the main sequence (Fig. [9]|. The bold parts of the curves indi- 
cate positive, outward fluxes that are associated with the coun- 
terclockwise meridional currents described in Fig. [6] 

Unlike in the massive star model, shear extracts angular mo- 
mentum everywhere because of the monotonous decrease of Q. 
with radius. Nevertheless, its contribution remains negligible. 
At the turn-off, the meridional circulation and angular veloci- 
ties have dropped in the outer parts of the radiative zone because 
of the braking of the surface layers and structural changes. This 
weakens Fmc while Fs remains almost unaltered. In the central 
regions, both Q and cJ^Q substantially increase leading to a con- 
comitant rise in F^c and Fs- The increase and sign change of 
Fmc between r - 0.05 Rq and r - 0.4 Rq is partly due to the 
off-center displacement of the nuclear energy production region 
when H is exhausted in the core. 

In summary, our diagnostic tools applied on two well-known 
cases allow us to give quantitative estimates on the main proper- 
ties of angular momentum transport already discussed in the lit- 



erature (see e.g., Meynet & Maeder|2000 Palacios et al.|2003] l. 
In low-mass stars, angular momentum transport is ensured by 
meridional currents that are primarily generated by the action of 
the applied torques resulting from the action of magnetic brak- 
ing. In the massive star, owing to the diff'erences in structure 
(density stratification, radius) and in surface velocities, the gra- 
dient of angular velocity becomes positive in the outer layers of 
the star and there shear turbulence takes over the transports of 
angular momentum. 



3.4. Thermal relaxation 

Let us now analyze how thermal relaxation is achieved, as de- 
scribed by the heat equation ( [T8| ), which rules the temperature 
perturbations. As explained previously, the meridional circula- 
tion is mainly driven by the extraction of angular momentum due 
to angular momentum losses at the surface (i.e. by the external 
torques) and to structural readjustments (such as the expansion 
of the star on the MS). This circulation advects specific entropy, 
which perturbs the thermal equilibrium and induces temperature 
fluctuations described by ^'2. These are linked to the differen- 
tial rotation through the baroclinic equation ( 



17 1, so that the heat 



equation is actually tightly coupled with the equation of angular 
momentum transport ( [TT] i. We examine how this works in our 
models. 



3.4.1. 20 Mr model 



Each term of Eq. ( 19 1 is represented in Fig. 10 (left panel). The 
plot is separated in two panels, the upper and lower ones show- 
ing the positive and negative parts of each profile respectively, on 
a logarithmic scale. The differential rotation adjusts itself so that 
the advection of entropy by the meridional circulation (TAdv) 
is almost exactly compensated by the thermal readjustment via 
the barotropic (Ts) and thermic (Trh) terms. The other terms, 
associated with the non-stationarity (Tns) and with the nuclear 
and gravitational energy generation (Tn-g). are several orders 
of magnitude smaller and thus can be neglected during the main 
sequence evolution. From this simplification, splitting the diver- 
gence of the radiative flux (third term in the right-hand side of 
Eq. 18 1 into barotropic and thermic (baroclinic) components thus 
appears to be somewhat artificial. 
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Fig. 10. Representation of the different components of the heat equation (Eq.[l8]l for the 20 Mq model when - 0.32 (left panel). 
The upper panel shows log|S| for S > while the lower panel shows log|S| for S < 0. Shown are the non-stationary (NS, dot 
long-dashed magenta line), barotropic (S, dotted blue line), baroclinic (Th long-dashed red line), advective (dashed black line) 
terms and the contribution due to the nuclear and gravitational energy variations (N-G, solid green lines). The central panel shows 
the components of the barotropic term; (solid line), Tgi (dotted line), 7^3 (long-dashed line), 7^4 (dotted-dashed hne). The 
right panel shows the thermic components: Ttii (solid), Ttiii (dashed line) and Trhi (dotted line). 
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Fig. 11. Two-dimensional reconstruction of the temperature fluc- 
tuations ST = r^2^'2(cos 0) in the 20 Mq model when X^. - 
0.32. The light (blue) parts are for negative fluctuations and the 
darker (red) are for positive ones. 



The temperature perturbations quickly relax to the 
asymptotic regime, which involves only the source (TAdv) and 
diffusion (Tg + Tjh) terms (see also Appendix B for description 



of these terms), and Eq. ( [T8] l gives: 

V ■ OfVr) -V-Fh^ -pTUrdrS 

and 



T2,S+T2,TU^pTCp^U2 



(31) 



(32) 



To explicitly disentangle the role played by each component 
in shaping the entropy transport, we split the barotropic term as 
follow: 



where 



7s 1 



.L 2 
'm3 



■2 \ 



2nGp 



—2 / 

Q dr 



and 



L 4 p„ — 2 / r 

M 9 p ' ~ 

. _L2 (e + ^gmv) -2^ (r\ 



tm = -p~^('pA2-6'¥2). 
M 3 p 



(33) 
(34) 
(35) 
(36) 
(37) 



7~si represents the "pure" barotropic term (which does not de- 
pend on 5, Q), while Tgi and Ts4 are the first and second baro- 
clinic corrections due to the differential rotation (cf. Eq. 17 i 
Finally, Tsi is associated with the energy production on the iso- 
bar 

Next, the thermic term is decomposed as follows: 



Trh - 7~Thl + 7~T| 



(38) 
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Fig. 12. Same as Fig. 10 for the 1.5 Mq model at the three evolutionary points indicated by X^. All quantities are shown in 



where 



and 



M\p3 3r 



_L2HtDi, 

-p Tt- 



(39) 



(40) 



The term Txhi comes from the diffusive laplacian acting on 
while '7Th2 is related to the heat flux carried by the horizontal 
turbulence (V ■ Fh). 

Middle and left panels of Fig. [10] display these components 
for the barotropic (Tg) and thermic (Txh) terms. We first note 
that the behaviour of Tg closely follows that of Tgi . This result 



is expected as the terms Tga and 7^4 are driven by the T- and 
//-fluctuations, which remain small in main sequence stars (see 
below). Moreover the termTga scales with (e+egi-av)/fm (i-C-, the 
ratio of energy production at a given radius with respect to the 
total energy enclosed in the sphere of this radius), which remains 
small as the (nuclear) energy production is mainly located in the 
convective core. It should also be noted that the relative strength 

of the two baroclinic corrections scales with Q drl^j^ which 

is always very small (between 10"'"* and 10"'" in the whole ra- 

—2 

diative region and mainly varies via the D. dependence). This 



explains why the term Ts2 does not appear in Fig. 10 (middle 



panel) and can be neglected during the main sequence evolution. 
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Fig. 13. The temperature fluctuation 6T - T^2P2(cos6) in the 1.5 Mg model at three different points on the main sequence, when 
the hydrogen mass fraction in the center = 0.675 (left panel), - 0.32 (middle panel) and Xc = (right panel). 



The thermic term (Txh) requires more attention as no specific 
component dominates in the whole radiative layer. Below the 
surface and above the convective core the term Txhi dominates 
while Txh is driven by "Fxha in the center of the radiative zone. 
As it will be discussed in the next section, this result strongly 
depends on the DhlKj ratio. 

Figure[TT]presents the two-dimensional reconstruction of the 
temperature perturbations in the meridional plane, i.e. 6T (r, 6) - 
^2(r)T(r)P2(cos6). These are directly connected to the be- 
haviour of Q-gradient: the temperature is higher at the poles 
(closer to the rotational axis) than at the equator where dr^l < 0, 
i.e. where the angular velocity decreases with increasing radius. 
Near the surface, the situation is reversed since S^Q is positive 
there. 



3.4.2. 1.5 Mo model 



stationary terms become dominant in the deep interior and must 
not be neglected anymore for the following evolutionary phase. 



Figure 12 is similar to Fig. 10 for the 1.5 Mq model, but con- 
trarily to the 20 model, 7"Adv is now of the same order of 
magnitude as and "Fxh. This is because magnetic braking 
forces a larger extraction of angular momentum. Let us note that 
the reversal of the meridional circulation shown in Figs. |6] and |7] 
is evidenced by the sign change of TAdv Also in contrast with 
the massive star, the nuclear and gravitational term, Tn-g, be- 
comes comparable to the advective one near the convective core, 
and should thus not be neglected. Finally, except near turn-off 
where it becomes important in the interior, the non-stationary 
term, Tns, is generally several orders of magnitudes smaller than 
the other terms during most of the main sequence. 



Equation (19 1 is thus left with four terms balancing each 



other. The thermic (baroclinic) and barotropic terms have similar 
amplitudes and opposite signs and their sum balances the sum of 
the advective term and the nuclear and gravitational heating: 



7~S -I- Tjh - Tadv + "^N-C 



(41) 



Let us note finally that this situation is not valid anymore at the 
turn-off' as seen in the third column of Fig. [12] There, the non- 



Figure 12 also displays the components of the barotropic and 
thermic term in the same way as in Fig. [10] In most of the radia- 
tive interior follows the barotropic component, Tgi , except 
below the surface where the contribution of the second baro- 
clinic correction becomes relevant. As in the 20 Mq model, 
Tsa is negligible during the main sequence. However in the last 
model (Xc - 0.0), this term dominates in the center due to the 
release of (a) nuclear energy by the hydrogen-burning shell and 
(b) gravitational energy by the contracting core. 

Contrary to the 20 Mq, the heat diffusion, Txhi (°^ DhlKj), 
plays no role in the thermic term, Txh, which is solely driven by 
the laplacian component, T'jhi . The 1 .5 Mq model has been com- 
puted with the Zahn (1992) prescription for the horizontal turbu- 
lent diffusion coefficient Dh leading to D/JKt - 10 - 10"^, 
whereas this ratio amounts to about 10 for the 20 Mq, where the 
Mathis et al.|(|2004 ) prescription was used (see § 3.5). At the cen- 
ter of the last shown model, the laplacian component presents a 
very complex behaviour as it is sensitive to the first, second and 
third radial derivatives of Q. Fig. [3] reveals that the profile of Q. 
presents some rapid variations in this region which explain this 
intricate feature. 

Figure[T3]displays the temperature fluctuations which essen- 
tially follow the gradient of Q as in the 20 Mq case. The main 
difference is due to the presence of a second circulation loop. 
The amplitude of the fluctuations remains globally of the same 
while the star is on the MS, and amount to about a few percents. 



3.4.3. Possible simplifications 

From the previous analysis of the 1.5 and 20 Mq models, during 
the quiescent main sequence evolution, the relaxation of *1'2 is 
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pTCp^U2. 



(42) 



Note that T^a is now expressed directly as a function of Q. where 



Eq. ( 17 1 has been used. On the other hand, we have verified that 
the A2-terms can be neglected compared to the ^2 ones during 
the main sequence evolution which allows simplifications in the 
writing of Txhi ■ The advection of entropy is then balanced by 
Tsi , 7^4 and the thermal diffusion. The simplification given by 
Eq. ( |42) i is no longer valid if structural readjustments occur, as 
it is the case when the core contracts at the end of the main se- 
quence^ 
If 



Zahn 



( 1992 I prescription for D/, is used. 



be furthermiore omitted. 



the term -p^ can 



n conclusion, the picture of the heat transport presented 
here differs radically from the classical Eddington-Sweet vision 
where the thermal imbalance induced by the centrifugal accel- 
eration is considered to be at the origin of the meridional cir- 
culation. Starting from an arbitrary rotation and thermal state, 
the advection of entropy leads, after thermal relaxation, to a new 
temperature state and thus to a new differential rotation profile 
(cf. EqfTTl). There then remains a circulation induced by the an- 
gular momentum losses and the structural adjustments, and the 
concomitant evolution of the rotation and thermal states. 



3.5. Transport oi chemicals 

Once a composition gradient builds up, due to gravitational set- 
tling (not accounted for in the present calculations) or nuclear 
burning, the meridional circulation creates a latitudinal pertur- 
bation of the molecular weight distribution, as described by 
Eq. (27 I. The sign of these /^-fluctuations depends on the ori- 
entation of the circulation, and we thus expect some differences 
between the two models. 



3.5.1. 20 Mq model 

In that model, the meridional circulation always transports mat- 
ter with low-ju down to the interior near the rotational axis, while 
high-ju matter is carried outward along the equator. As a result, 
a horizontal ju-gradient appears near the convective core (see 
Fig. [14]), which is directed from pole to equator, i.e. opposite to 
that of the temperature fluctuations. These yu-fluctuations remain 
small (less than 1%), which justifies the perturbative approach. 

We already mentioned in § |2.3.4| that the strong turbulent 
transport in the horizontal direction interferes with the merid- 
ional circulation, rendering the vertical transport of diffusive na- 
ture, characterized by an effective diffusivity Deff (C haboyer &| 
[Zahn 1992). This description makes the comparison between the 
two processes that contribute to the transport of chemical el- 
ements (vertical shear and combination of advection and hori- 
zontal shear) easier since it suffices to compare their respective 
diffusivities. 



This is done in Fig. 15 where we see that in most of the stel- 
lar interior the vertical shear turbulence largely dominates over 




-8.21e-04 



4.10e-04 



Fig. 14. Two-dimensional reconstruction of the //-perturbations, 
612 - JIA2P2 (cos 6\ in the 20 Mq star when - 0.32. Hatched 
regions delineate convection zones. 




2 4 6 8 
r (Ro) 

Fig. 15. Profiles of the thermal diffusivity (Kj dotted-dashed 
cyan line) and that associated to the meridional circulation (Dgff 
short-dashed red line) and with the horizontal (D/, long-dashed 
black line) and vertical (Dj, dotted blue line) turbulence in the 
20 Mq model when = 0.32. 



the combination of advection and horizontal shear (D,, » Oeif), 
except in the vicinity of the convective core, where the strong 
/(-stratification inhibits the shear instability, and in the region 
where the gradient of Q changes sign (and where D,, reduces 
to the molecular diffusivity). 



This figure can be compared with Fig. 6 of Meynet & Maeder 
( |2000| l, where the same profiles are presented for a 20 Mq star 
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atZ- 0.02, with initial equatorial velocity Vi^i = 300 km s ' . In 
their simulation, Meynet & Maeder use the prescription given by 
Eq. ( 30 1 for the horizontal turbulent viscosity. While we obtain 
very similar profiles for the vertical shear diffusion coefficient, 
Dv, and for the thermal diffusivity, Kt, the diffusion coefficients 
associated with horizontal turbulence, Di,, and meridional circu- 
lation, ZJefF, differ by several orders of magnitude. The Mathis| 



et al. ( 2004 1 prescription results in larger horizontal turbulent 
viscosity, which also translates into lower ZJgff. This very same 
result was firstly obtained by |Maeder| ( |2003| l, where using ener- 
getic considerations, he derived a new expression for D/, that is 



very similar to Eq. 29 (see also Mathis et al. 2007). It should be 
noted that the anisotropic turbulence assumption (D/, » D,,) is 
always satisfied. 



structural changes external torques 



meridional circulation • 



baroclinic temperature field 

Baroclinic 
torque 



differential rotation ► shear turbulence 



3.5.2. 1.5 Mo model 

Figure [16] shows that the diff'usion coefficients in the 1.5 Mo star 
are much smaller than in the massive star as a result of lower an- 
gular velocity and degree of differential rotation. The diffusion 
coefficients Deft and D,, are much closer in this model, partly 
due to the adoption of Zahnjs ( |1992 1 prescription for D/, (see 
Eqs. (|29]) and ([30]l and Mathis et al. 2004 for more details). 
However, the decreasing efficiency of the meridional transport 
as the star evolves on the main sequence results in the shear tur- 
bulence dominating the distribution of chemicals (see Palacios] 
et al. 112003 I. The inhibiting action of the composition gradient 
is visible near the convective core where D,, drops to its micro- 
scopic value, indicating that turbulence is suppressed. It should 
also be noted that in this model the anisotropic turbulence as- 
sumption (D/, » Dy) is also always satisfied. 

The sign of the /i-perturbation depends on the orientation 
of the circulation near the convective core. In the lower mass 
model, differences arise from the development of a secondary 
loop, as seen at = 0.675 and X^. = 0. 

Although the meridional circulation and Di, are both smaller 
in the low-mass star, the ratio U2/D1,, which governs the rela- 
tive strength of vertical advection by meridional circulation and 
of horizontal diffusion (see Eq. 27 1, is larger in the 1 .5Mo com- 
pared to the 20 Mq. This results in larger yu-imbalance in the 
lower-mass model. 



4. Conclusion 

The paper presents a complete set of diagnostic tools that pro- 
vide a comprehensive and coherent understanding of the secu- 
lar hydrodynamical transport processes operating in the radia- 
tive zones of rotating stars. The framework is that of "type I 
rotational mixing", where angular momentum and nuclides are 
transported by the same mechanisms, namely large-scale merid- 
ional circulation and shear-induced turbulence. 

In order to validate our new approach, the first analysis pre- 
sented here is performed on two well-studied cases of a low- 
mass (1.5 Mq) (see Palacios et al.|2003) and a massive (20 Mo) 
star ( |Meynet & Maeder|2000 1. It allows us to disentangle what 
are the main processes at the origin of meridional circulation and 
driving angular momentum, heat and chemicals transports in the 
stellar interior and to confirm results previously established in 
the literature. In particular, the angular momentum loss, either 
by radiation-driven wind (the massive star) or by magnetic brak- 
ing (the low-mass star) combined with structural readjustments 
are the sources responsible for the generation of meridional cir- 
culation (Fig. 18 1. The direction of the flow is such that it trans- 
ports angular momentum outward to the surface, except in the 



Fig. 18. Rotational mixing of type I in the radiative zone of a ro- 
tating star, where the transport of angular momentum is achieved 
by meridional circulation and turbulent diffusion. 



deep interior of the 1 .5 Mq star where it can reverse, producing 
a steeper angular velocity gradient. 

This circulation advects heat, and thus generates latitude de- 
pendent temperature fluctuations; these tend to be damped out 
by radiative diffusion, until they establish a subtle balance be- 
tween advection and diffusion, that allows them to induce differ- 
ential rotation through the baroclinic torque (Eq. [TTJi. This dif- 
ferential rotation is shear-unstable, and generates turbulence that 
paiticipates to the transport of angular momentum, thus clos- 
ing the loop (Fig.frsj). This picture thus differs drastically from 
the classical Eddington-Sweet circulation, where the circulation 
is deemed to originate from the thermal imbalance due to the 
centrifugal force. If turbulent transport and extraction of angular 
momentum were absent, the circulation would die altogether, as 



was pointed out by Busse ( 1982 1. 



Thanks to the diagnostic tools we have developed, the 
present investigation helped to clarify a number of other points. 
We confirm that it is mainly the meridional circulation that trans- 



ports angular momentum in the low-mass star (see also^Palacios 
|et al.|20(j3] l; it is thus crucial to describe that transport as an ad- 
vective process, as in the approach developed by Zahn, Maeder, 
and collaborators. In the massive star, where the redistribution 
of angular momentum due to the star's expansion plays a much 
more important role than its extraction by the radiative wind, ad- 
vection and diffusion make comparable contributions, though in 
different regions: the advective transport dominates in the inner 
layers of the star, whereas the diffusive transport is the prime ac- 
tor in the upper part and, moreover, works there in the opposite 
direction. 

For the first time, we confirm through detailed computation 
and separation of the different components, that the advection 
of entropy by meridional circulation is almost exactly balanced 
by the thermal relaxation (barotropic and thermic terms) dur- 
ing most of the main sequence. This was generally assumed 
prior to the [Maeder & Zahn (1998) work integrating the time- 
dependency of the temperature fluctuations, but not actually ver- 
ified through the computation. Grounding our detailed analy- 
sis on advanced graphical tools, we are also able to propose a 
simplified expression for the thermal relaxation on the main se- 
quence. 

On the other hand, the transport of chemicals is shown to 
be dominated by the shear-induced turbulence while the star 
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Fig. 16. Same as Fig. 15 for the 1.5 Mq model when = 0.675, = 0.32 and X^ = 0.0. 
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Fig. 17. Two-dimensional reconstruction of the yU-perturbations, 5^ — IJ.IS.2P2 (cos 6) in the 1.5 M© when — 0.675 (/e/f panel), 
= 0.32 {middle panel) and = 0.0 (Wg/if panel). Hatched regions delineate convection zones. 



evolves on the main sequence. This is particularly clear for a 
massive star such as the 20 Mq star 

In the near future, we plan to apply the same tools to Type II 
rotational mixing, where we shall include the transport of angu- 
lar momentum by internal gravity waves and magnetic stresses 
and to advanced phases of stellar evolution. 
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Appendix A: Meridional circulation and angular 
momentum transport 

Let us consider the Eq. ( [TT[ ) for the angular momentum transport: 
(r^Q) ^j^d, (pr^Qf/j) + -^5, (pv.-r^^.Q) . (A. 1) 



Integrating over a spherical shell, we get: 

J r=i-\ J r-r\ 

nr=r2 



(A.2) 



ri and r2 being respectively its inner and outer radius. Then, in- 
troducing the elementary mass element dm = 4npr^dr, we ob- 
tain: 



1 r-m=m{r2) J 1 r 



(A.3) 



We apply the following identity 



d 
df 



■mit) 

f(m, t)dm 

m,(t} 

(dmiit) 



■'«2(') ^ 

df 



[f{m, t)] dm 



mi (I) 



\ dt 



dmi(f) 

f(m2,t) finiut) 

dt 



In our Lagrangian description, the mass is conserved, 



(A.4) 

dm2(0 _ 



dm, (I) 
di 

1 d 

47T dt 



= 0, so 



Urn 



•m=m(r2) 



r'^Qdm 



m=m(ri ) 



(A.5) 



To obtain the value of the meridional circulation at m(r), we set 
ri = and r2 - r: 



U2 = [r(m) -pv,/^,^] 

^.Ar^ t J 



where 



nm(r) 

Jo ' 



r'^ndm' 



(A.6) 



(A.7) 



Appendix B: Components T2 of the heat equation 



These terms intervene in the heat relaxation equation ( [T8[ l and 
their derivation may be found in Zahn (1992) or Mathis & Zahn 
(2004). They include : 
the barotropic term: 



1 - 



Q. 



2p„; 



2;7rGp 3 p 

i^A2-6^2), 



^P'"t A AVr/^ ^ ^g"^) 
iipA2 - 6^2) 



n dr\ — 



3 p 

the thermic term: 
P. 



T2 



Th 



(B.l) 



(B.2) 



the term associated with local energy sources: 

(e + igrav) 



[^2 ir) + (Ut - feS + 
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where 

^2 W = HTdr'¥2 -{1-6 +XTm - W +Xu,)^2- (B.3) 

A^^ = ig6/Hp) (Vad - V) is the square of the thermal part of the 
Briint-Vaisala frequency where V = (91n r/(91nf is the radiative 
temperature gradient and Vad is the adiabatic one. L, M, Cp have 
their usual meaning and T is the horizontal average of the tem- 
perature. We have also introduced the pressure and temperature 
scale-heights Hp - |dr/dlnf | and Hj - |dr/dlnr|, the thermal 
diffusivity Kj = xIpCp = 16(tT^ /3Kp^Cp the horizontal eddy- 
diffusivity Dh and = e/(^? + £gnn), with ? and igrav being the 
mean nuclear and gravitational energy release rates respectively. 
cr and k are the Stefan constant and the Rosseland mean opacity. 

- (51ne/c?lnyu)pj- and;^'^ - (dlnx /din ij.)p-j. are the logarith- 
mic derivatives of e and of the radiative conductivity x with re- 
spect to ju, and £t = (dlne/dlnT)pjj stadxT = idlnx/dlnT)p^ 
are the derivatives of these same quantities with respect to the 
temperature T. Finally = L/M and pm = 3g(r)/4nrG are the 
horizontal average of the energy production rate and the mean 
density inside the considered level surface. 



